Quantum Hall effect and Landau level crossing of Dirac fermions in trilayer graphene 
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We investigate electronic transport in high mobility (> 100,000 cm 2 /V-s) trilayer graphene (TLG) 
devices on hexagonal boron nitride, which enables the observation of Shubnikov-de Haas oscillations 
and an unconventional quantum Hall effect. The massless and massive characters of the TLG 
subbands lead to a set of Landau level crossings, whose magnetic field and filling factor coordinates 
enable the direct determination of the Slonczewski-Weiss-McClure (SWMcC) parameters used to 
describe the peculiar electronic structure of trilayer graphene. Moreover, at high magnetic fields, 
the degenerate crossing points split into manifolds indicating the existence of broken-symmetry 
quantum Hall states. 



Bernal or ABA stacked TLG (Fig. lb) is an intriguing 
material to study Dirac physics and quantum Hall effect 
(QHE) because of its unique band structure which, in the 
simplest approximation, consists of massless single-layer 
graphene (SLG) and massive bilayer graphene (BLG) 
subbands at low energy (Fig. lc)[l|-|4|. The energies of 
the Landau levels (LLs) for massless charge carriers de- 
pend on the square root of the magnetic field \[B 
while for massive charge carriers they depend linearly on 
B @, §, Eq|- Therefore, the LLs from these two differ- 
ent subbands in TLG should cross at some finite fields, 
resulting in accidental LL degeneracies at the crossing 
points. 

However, one of the major challenges so far to ob- 
serve QHE in TLG has been its low mobility on SiC>2 
substrates [ll|,[l2|. To overcome this problem, we use high 
quality hexagonal boron nitride (hBN) single crystals [13j| 
as local substrates, which have been shown to reduce 
carrier scattering in graphene devices [14j. Substrate sup- 
ported devices also allows us to reach higher carrier den- 
sity than suspended samples [HI, which is necessary for 
the observation of the LL crossings. 

Figure la shows an atomic force microscope image 
of a Hall bar shaped TLG device on hBN. Our fabri- 
cation process consists of mechanically exfoliating hBN 
and graphene flakes on different supports, and a flip chip 
bonding step to align them on top of each other (see 
SI for details). The graphene flakes are then patterned 
into a Hall bar geometry and contacted by electron beam 
lithography. The device is then annealed in forming gas 
to remove residue and cooled down in a He-3 cryostat. 

In order to further reduce disorder and increase 
the mobility, we perform current annealing at low 
temperature [16|. Figures le and If show the resistiv- 
ity and conductivity of a TLG device at zero magnetic 
field after current annealing. The resistivity at the Dirac 
peak exhibits a strong temperature dependence, which in 
SLG is a strong indication of high device 

quality mEIt 

In addition, we also observe a double-peak structure at 
low temperatures (Fig. le). This double-peak structure 
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FIG. 1: a) False color atomic force microscopy image of a TLG 
Hall bar device on hBN. b) Bernal stacked TLG atomic lat- 
tice. The SWMcC hopping parameters, 7*, are shown by pur- 
ple dashed lines connecting the corresponding hopping sites. 
In addition to ji , the SWMcC parameters also include the on- 
site energy difference, 5, between A and B sublattices (blue 
and red lattices), c) Band structure of TLG at low energy, 
which takes into account only the nearest neighbour intra- and 
inter-layer hopping parameters 70 and 71 . d) Band structure 
of TLG within a full parameter model, with the parameters 
calculated from the SdH oscillations in Fig. 2b. e) Resistiv- 
ity as a function of density and temperature for TLG. The 
double peak structure starts to emerge as temperature de- 
creases below 10 K. f) Conductivity as a function of density 
and temperature. The field-effect mobility at 300 mK reaches 
-110,000 cm 2 /V-s and decreases to -65,000 cm 2 /V-s at 40 K 



is likely due to the band overlap which occurs in TLG 
when all SWMcC parameters are included in the tight- 
binding calculation of its band structure, as we show 
below. The field effect mobility of this device reaches 
110,000 cm 2 /V-s at 300 mK at densities as high as 
6 x 10 11 cm -2 . This mobility value is two orders of mag- 
nitude higher than previously reported values for sup- 
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FIG. 2: a) p xx as a function of inverse magnetic field at 300 mK. The numbers inside the figure indicate the filling factors at the 
SdH oscillation minima. The highlighted bands show the regions of 8-fold degeneracy, which provide evidence for LL crossings 
of the SLG- and BLG-like subbands. For B > 4 T, the SdH minima are separated by Ais = 1 or 2, indicating the splitting 
of LLs. b) Color map of p xx versus n and B at 300 mK. The diagonal lines correspond to constant filling factor lines. The 
beating pattern, most visible at negative densities, is a consequence of LL crossings. The white central region corresponds to an 
insulating state at zero density (see SI), c) Calculated LL energy spectrum in TLG for the SWMcC parameters obtained from 
b. The red dashed and black lines are LLs at K and K' points respectively. The roughly V^B-like and linear B-like dispersion 
from the SLG- and BLG-like subbands is evident. Each line corresponds to a spin degenerate LL. d) Calculated density of 
states as a function of density and B from the LL spectrum in c. Apart from the LL splitting, the location of the LL crossings 
agrees very well with the experimental data in b. e) p xx and cr xy as a function of filling factor for B — 7.3, 7.8, and 8.3 T. The 
highlighted orange region shows the appearance of the LL crossing at v — — 23 while the green highlighted region shows the 
LL splitting occurring at v — —20. 



ported TLG HiI, llil a nd comparable to suspended SLG- 
TLG samplesHa ll8| . The low disorder and high mo- 
bility enable us to probe LL crossings of Dirac fermions 
through the measurement of Shubnikov-de Hass (SdH) 
oscillations. 

Figure 2a shows longitudinal resistivity p func- 
tion of 1/5, for a carrier density n = —4.4 x 10 12 cm -2 . 
A pattern of SdH oscillations is clearly visible, albeit 
with different visibility and features depending on the 
B range. At low B (below ~ 1 T), there are a number 
of oscillations characterized by broad minima separated 
by relatively narrower maxima. Beyond ~ 1 T, the min- 



ima become sets of narrower oscillations, and a clear pat- 
tern emerges: each minimum in the oscillations indicates 
a completely filled LL with corresponding filling factor 
v = hn/eB, where h is Planck's constant, and e is the 
electron charge. Within a single particle picture, each 
LL is 4- fold degenerate, the degeneracy originating from 
the valley (K and K') and spin (up and down) degrees of 
freedoms in both the SLG-like and BLG-like subbands. 
When LLs from these two subbands cross at a given 5, 
the coexistence of two 4-fold degenerate LLs increases 
the degeneracy to 8-fold. This 8-fold degeneracy is high- 
lighted by the green bands in Fig. 2a, where v changes 
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by 8 from minimum to minimum instead of by 4. For 
B > 4 T, the splitting of the LLs results in v chang- 
ing by either 1 or 2, as the different broken- symmetry 
quantum Hall states are occupied. 

A more complete understanding of the TLG LL energy 
spectrum is obtained by plotting p xx as a function of n 
and B as shown in Figure 2b. The resulting fan dia- 
gram lines correspond to the SdH oscillations mentioned 
above, while the white central region corresponds to an 
insulating behavior at v = (see SI for details). The 
abovementioned crossings of SLG-like and BLG-like LLs 
manifest themselves as a beating pattern in the SdH oscil- 
lations, with a greater number of them and more visible 
on the hole side (n < 0). This electron- hole asymme- 
try results from the TLG band structure, as we show 
below. In addition, the LL splittings appear as finer 
split lines in the SdH oscillations. For each LL cross- 
ing, there is an enhancement of p xx due to the enhanced 
density of states [HI [lo|, and each crossing point can be 
uniquely identified by B and v. For instance, at B ~ 3 T 
and n ~ —4 x 10 12 cm -2 , the filling factors associated 
with the minima in the corresponding SdH oscillations 
change from 50 to 58 indicating that the crossing occurs 
at v = 54. 

The positions of the crossings in B and v space depend 
sensitively on the TLG band structure, and therefore en- 
able a direct electronic transport determination of the rel- 
evant SWMcC parameters for TLG. These parameters, 
proposed to explain the band structure of graphite (21). 
describe the different intra- and inter-layer hopping terms 
in the different graphene sheets (Fig. lb). We note that 
TLG is the fewest layer graphene system whose descrip- 
tion includes all the SWMcC parameters. The simplest 
TLG model, in which only the nearest intra- and inter- 
layer couplings (70 and 71) are considered (the ones typ- 
ically used to describe SLG and BLG), results in sym- 
metric electron and hole bands (Fig. lc) and therefore is 
clearly insufficient to explain the experimental data. We 
therefore use all the relevant SWMcC parameters to nu- 
merically calculate the LL energy spectrum (Fig. 2c) and 
density of states as a function of B (Fig. 2d), and perform 
a minimization procedure to fit the experimental data in 
Fig. 2b. In order to lower the number of parameters, we 
take 70 = 3.1 eV, 71 = 0.39 eV and 73 = 0.315 eV (see 
SI), and we obtain from our fit the following values of the 
SWMcC parameters; 72 = -0.028(4) eV, 74 = 0.041(10) 
eV, 75 = 0.05(2) eV, and S = 0.046(10) eV. The defini- 
tions of the ji can be found in Fig. lb and S is the on- 
site energy difference between the two-inequivalent car- 
bon sublattices residing in the same graphene layer. The 
values of the SWMcC parameters obtained are similar 
to previously reported values for graphite [2l| and, apart 
from the broken-symmetry states (see discussion below), 
our data agree very well with the LLs corresponding to 
Bernal stacked TLG, and not to rhombohedral stacked 
TLG [22]. These parameters result in the overall electron- 
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FIG. 3: a) p xx as a function of density and B at 300 mK 
showing a manifold of LL crossing points. The high p X x re- 
gions correspond to enhanced degeneracy due to LL crossings. 
Five crossing points are clearly visible and the sixth point is 
starting to appear in the lower-right corner. White dashed 
lines are guides to the eye for each v labeled on the edges, 
b) Schematic splitting and crossing of LLs yielding the mani- 
fold of crossings shown in a. Red and blue lines represent the 
split LL spectrum for the broken-symmetry QH states of the 
N = -ILL from the SLG-like subband and the N = -5 LL 
from the BLG-like subband, respectively. The degeneracies 
for each level are g — 1 for thin lines and g — 2 for the thick 
line. The highlighted green area corresponds to the region 
observed in the data in a. The numbers inside each region 
show the corresponding filling factors 



hole asymmetric band structure shown in Fig. Id, with 
small band gaps E gy s ~ 7 meV and E gy B ~ 14 meV, for 
the SLG- and BLG-like subbands, and a band overlap 
E Q ~ 14 meV. 

The LLs in TLG are not truly 4-fold degenerate even 
in a single particle picture, owing to the finite value of 72, 
75, and J, which break valley degeneracy (see Fig. 2c), 
in addition to the Zeeman interaction which breaks spin 
degeneracy. Our data at high B (Fig. 2a and 2b) show 
that the splitting of 4-fold degenerate LLs is observed up 
to filling factors as high as v = 46. While single par- 
ticle effects may partly explain these broken-symmetry 
QH states (e.g. from the width of the LLs crossings, we 
estimate the disorder broadening of the LLs to be ~1 
mV which is similar to the the Zeeman splitting at ~8 
T), it is likely that electron-electron (e-e) interactions 
play a significant role too, as it is the case in SLG and 
BLG[14|, l23l426j . For example, the insulating behavior 
we observe at v = 0, cannot be explained by single par- 
ticle effects, given the band overlap between the SLG- 
and BLG-like subbands, and the single particle LL en- 
ergy spectrum shown in Fig. 2c. However, a more de- 
tailed study including measurements of the gap energies 
and measurements in tilted magnetic fields, beyond the 
scope of this paper, is necessary to investigate the precise 
role of e-e interactions in TLG. Figure 2e shows example 
traces where the different behavior of LL crossings and 
LL splitting can be seen. 

At high £>, the LL crossing points should become cross- 
ing manifolds due to the crossing between the split SLG- 
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FIG. 4: cr xy and p xx as a function of density at B = 9 T 
and T — 300 mK, and before the last current annealing step. 
The dashed lines indicate the expected QH plateaus based 
on the simplest TLG model approximation. The dotted lines 
indicate the extra QH plateaus based on the full band struc- 
ture determined from Fig. 2c. (Inset) Calculated Density of 
states using full SWMcC parameter model. The blue line is 
calculated using higher disorder broadening than the red line. 



and BLG-like LLs. One such example is shown in Fig. 3a. 
From the LL energy spectrum shown in Fig 2c, the mani- 
fold corresponds to the crossing between the N = — 1 LL 
of the SLG-like subband, LLg 1 , and the TV = — 5 LL of 
the BLG-like subband, LLg 5 . In order to reproduce the 
observed degeneracies at the crossings, the 4- fold LLg 1 
has to completely split into four singly-degenerate LLs 
while the 4- fold LLg 5 splits into 3 LLs: two singly degen- 
erate LLs and one doubly degenerate LL. Figure 3b shows 
schematically the full 12-point manifold, of which only 6 
crossing points are visible in our density and B range. 
We have found that this splitting scheme is the only one 
that yields the correct result for both the degeneracies at 
the crossings and the filling factors at which they occur. 
The observation of the full 4- fold splitting of the LLg 1 
in TLG, although expected, is remarkable since previous 
transport studies of the N = 1 LL in SLG had reported 
only the breaking of some of the degeneracies 23, l27| . 
and the full 4-fold splitting has only been seen in re- 
cent STM experiments HH. The 1-2-1 splitting of LLs 
from the BLG-like subband, however, is more anomalous. 
Naively, one would expect the splitting to be either 2 fold 
or 4 fold, depending on whether one of the two degrees 
of freedom (valley or spin) is split or both are 25[ [26| . 
However, we note that this 1-2-1 splitting may also be 
present in a recent study of BLG on hBN in the inter- 
mediate ^-regime [li). and may possibly indicate a richer 
phase diagram based on SU(4) rather than SU(2)xSU(2) 
symmetry breaking. A detailed study of the crossing be- 
tween spin/ valley polarized LLs of massless and massive 
Dirac Fermions, together with the aforementioned possi- 
ble role of e-e interactions, could potentially lead to some 
intriguing phenomena such as phase transitions in quan- 
tum Hall ferromagnets ill [29| . 



Although the splitting of the LLs at high B provides 
insight into broken symmetries in TLG in the QH regime, 
it also masks out the QH plateaus expected within the 
simplest single particle model for TLG. The sequence of 
plateaus arising from such simple models has proven a 
useful tool in identifying SLG and BLG0, @, Ejj. For 
completeness, Figure 4 shows p xx and <j xy at B = 9 T 
before current annealing, i.e. in the presence of increased 
disorder which prevents the observation of LL splitting. 
In the simplest model, the QHE plateaus are expected at 
c^xy = ±4(N + 1/2 + l)e 2 /h for iV = 0,1,... where the 
12-fold zero energy LL results from the 4- fold and 8-fold 
zero energy LLs of the SLG- and BLG-like subbands, 
respectively [Sol O ur observations agree with this 

simple prediction for \v\ > 10 (with observed plateaus 
at ±10, ±14, ±18e 2 /fo), but we observe in addition ex- 
tra plateaus for v = ±2 and ±4 as well as the absence 
of a plateau at v = ±6. This unconventional QHE can 
be explained within the band model calculated using the 
SWMcC parameters obtained from Fig. 2a-c. In such 
model, the non-zero values of 72, 75, and S lift the degen- 
eracy of the "zero-energy" LLs of the SLG- and BLG-like 
subbands (Fig. 2c). In addition, the 4- fold degenerate 
N = LL of the SLG-like subband splits into two 2-fold 
degenerate valley polarized LLs and the 8-fold degener- 
ate (spin, valley and N=0,1 LLs) zero energy LLs of the 
BLG-like subband splits into two 4-fold degenerate LLs 
(the splitting between N = and N = 1 LLs remains rel- 
atively small compared to the valley splitting). We note 
that the Zeeman splitting is at least an order of magni- 
tude smaller than other types of splitting even at 9 T 
which is the reason why LLs remain spin degenerate in 
this non-interacting model. 

The inset to Fig. 4 shows the calculated density of 
states as a function of energy at 9 T. The zero density is 
located between two nearly degenerate LLs, each with 2- 
fold degeneracy which explains the observed plateaus at 
v = ±2. The absence of a plateau at v = is likely due to 
disorder, which smears out the small energy gap between 
these two LLs. The plateaus at v = ±4 stem from the 
next 2-fold degenerate LLs. However, these plateaus are 
not yet completely developed at 9 T, especially the one 
at v — — 4 (<j xy = 4e 2 /h) which coincides with the small 
energy gap between this LL and the next one. Finally, 
the absence of a plateau at v = ±6 (<r xy = —6e 2 /h) is due 
to the crossing between a 2-fold and a 4-fold degenerate 
LL. The degeneracy at the crossing becomes 6-fold and 
causes the position of the plateau to step from v = 4 to 
v = 10 (the non-developed v = 4 plateau does not reach 
its exact value at <7 xy = —Ae 2 /h). Unlike SLG and BLG 
in which the sequence of the plateaus are the same for all 
£>, the observed plateaus in TLG depend on B because 
of the LL crossing. 
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